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A Bose-Einstein condensate with attractive interaction can 
be metastable if it is spatially confined and if the number of 
condensate bosons No is below a certain critical value N c . By 
applying a variational method and the instanton techinique 
to the Gross-Pitaevskii energy functional, we find analyti- 
cally the frequency of the collective excitation and the rate of 
macroscopic quantum tunneling (MQT). We show that near 
the critical point the tunneling exponent vanishes according 
to (1 — No/N c )i and that MQT can be a dominant decay 
mechanism of the condensate for No very close to -/V c . 
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An explosively increasing amount of research is being 
carried out on the phenomenon of Bose-Einstein conden- 
sation (BEC) in gases of 87 Rb §, 7 Li §, and 23 Na |] 
atoms. Unlike the other two species, a uniform system of 
7 Li atoms is usually believed not to form a stable BEC 
state H because the s-wave scattering length a is negative 
and the attractive interaction between the atoms causes 
the condensate to collapse upon itself. When atoms are 
spatially confined, however, they acquire zero-point en- 
ergies which, under certain conditions, counterbalance 
the attractive interaction, thereby allowing a metastable 
condensate to form. 

As the number of condensate bosons increases, the at- 
tractive interaction becomes strong and the energy bar- 
rier that prevents the condensate from collasing becomes 
accordingly low. Given a confining potential which de- 
termines the zero-point energies, there exists a criti- 
cal number N c of condensate bosons at which the en- 
ergy barrier vanishes. When the number of condensate 
bosons Nq is slightly below N c , the energy barrier will 
be so low that the condensate might undergo macro- 
scopic quantum tunneling (MQT) to a dense state. Ka- 
gan et al. j^] estimated the overlap integral between the 
metastable condensate and the dense state to be pro- 
portional to exp(— |iVo In jt), where Iq is the amplitude 
of zero-point oscillations of the trap, and L* ~ 3|a| AT . 
Shuryak || estimated the MQT rate to be proportional 
to exp [-0.57(JV C - No)]. Stoof § wrote down a WKB 
formula for the MQT rate but did not explicitly evaluate 
it near the critical point. 

In this Letter we use a variational method and the 
instanton technique to show that the tunneling expo- 
nent vanishes faster than found in Refs. |^,^|, namely 



as (1 — Nq/N c )*, as Nq approaches N c . By comparing 
the MQT rate with other possible decay mechanisms, we 
argue that MQT can be a dominant decay mechanism of 
the condensate near the critical point at zero tempera- 
ture, contrary to the conclusions of Refs. [pLpI- Since the 
formulas we obtain contain no fitting parameters, they 
can be used as stringent tests of existence of BEC. 

At sufficiently low temperatures, a condensate of 
weakly interacting bosons is described by a single wave 
function 'J'(r), and the interaction between them is de- 
scribed by the s-wave scattering length a. The wave func- 
tion is determined so as to minimize the Gross-Pitaevskii 
energy functional Jq] 
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where m is the atomic mass, V(r) is a confining potential, 
and 2-kTi 2 a\^> (r)| 2 / 'm is the local mean- field interaction 
energy per particle. The wave function if? is normalized 
so that J |\&| 2 dr is equal to the number of condensed 
bosons Nq. We assume an axially symmetric confining 
potential: V(r) = m(uj\x 2 + uj\y 2 + lo 2 z 2 )/2, where uj± 
and are respectively the frequency of the radial con- 
fining potential and that of the axial one. 

To evaluate we assume a Gaussian trial wave 

function [01 
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where d± and d\\ are variational parameters. Substituting 
Eq. (||) into Eq. (||), we obtain 



= 2\-i (ri + rj 2 ) + \i (rjj + 
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where r± 



d±y/ mujj_/h 



d± I d±o and r 



d\\yJmu)\\/Ti = dii/diio are the radial and axial widths of 
the condensate normalized by their noninteracting val- 
ues; A = u)\\/u>± is the asymmetry parameter of the con- 
fining potential, and 
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is the dimensionless strength of interaction 

For a metastable condensate to exist, the function 
f(r± , r» ) must have a local minimum which is determined 
from conditions 



df_ 

dr± 
df_ = 

dr\\ 



= 
-> 



A37 



2(1-4)' 
7^|| 

2A§(l-rjf)" 



(5) 
(6) 



The mestability of the condensate is determined from the 
curvatures of f(r±,ru) which are calculated to be 



d 2 f 1 d 2 f 

-4 = 16A-*>0, 



d 



2Af(3 + r7 4 ) > 0, (7) 
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The condition for a metastable condensate to exist is that 
the curvature of f(r±,r\\) at the local minimum is posi- 
tive in all directions, that is, 



dr\ dr 2 



d 2 f 
dridr 



> 0. 
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It is clear from Eqs. (|t|) and (||) that the condition (^) is 
always satisfied at the limit of weak interaction 7 — ► 0, 
but that it is violated at the limit of strong interaction 
7 — ► 00. In between there must be a critical value of 7 
such that the left-hand side of (ph is zero. It is given by 



7c 
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For 7 > 7c the local minimum becomes a saddle point 
and the condensate becomes unstable, collapsing into a 
dense state. Substituting Eq. (1Q) into Eqs. (0) gives 



A*(l 



V 2 1c( 1 + 3 1c)' 
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where rii c denotes the value of r\\ at the critical point. 
Substituting Eqs. @ and (0) into Eq. (§) gives 



A 2 = 



4ic(l + 3rf c ) 2 



(1- r| c )3(3 + 54)- 

We may use Eq. ( |l2| ) to simplify Eq. (|l(J) somewhat 
= 4r|| c (l + 3rp 
7c A*(3 + 5r| c )' 



(13) 



For a given asymmetry parameter A = u>u/u>±, a real 
positive root of Eq. ( |l2| ) gives f|| c . Substituting this into 
Eqs. ( pT| ) and (|l^) gives r± c and rf c , respectively. 



For an isotropic case (A = 1), we obtain r± c = rn c = 
5-3 ~ 0.669 and 7 C — 1.07 in agreement with the results 
of Refs. 0,0] • Taking experimental data from the sec- 
ond reference in 0), where a = — 14.5A, d±o ~ 3.08/im 
and A ~ 0.867, we obtain from Eq. (|) that A^ c ~ 1460. 
This is 17% greater than the more precise value of 1250 
which is obtained by numerically solving the nonlinear 
Schrodinger equation [jll]]. In evaluating the MQT rate, 
we will use the latter value for the critical number of 
condensate bosons. 

Figure [l] shows r_L C , r|| c , and 7 C as a function of the 
asymmetry parameter A. We see that the maximum value 
of 7 C can be attained for the case of an isotropic potential 
p"2| . Also plotted is the ratio d± c /d\\ c of the width of 
the condensate along the radial direction to that along 
the axial one. We note that the ratio remains relatively 
constant for A < 1 but it grows rapidly for A > 1. In what 
follows we will focus on the case of an isotropic confining 
potential, and drop the subsreripts _L and || . The function 
/ can then be written as f(r) = 3r 



Collective excitation of the condensate 



The condensate undergoes density oscillations around 
the local minimum. Associated with this collective mo- 
tion is a kinetic energy T which we may write as 

T = 3CN md 2 = ZQNomdlr 2 = 3^—^r 2 , (14) 

where £ is a constant to be determined below. 

The dynamics of the collective excitation is determined 
by the kinetic-energy term (^) and the quadratic part 
of the potential energy which is obtained by expanding 
/( r ) — /( rlnm ) m powers of r — r mm : 



^/V!) (r rm in ) 2 ) 



(15) 



where r mm is determined from /'(r mm ) = 0. For 7 
slightly below 7 C , r mm is close to r c = 5~«; that is, 
r mm = r c + 5r, where Sr is given to leading order in 

7c - 7 by 



S r ~ _f_( 7c _ 7 )j 



(16) 



(12) ^Frorn Eqs. (|) and @), /"(r™ 1 ") is calculated to be 



/"(ri m ) =s 120 . 

T I r. 



(17) 



^From Eqs. (|l4|)-([l7|) we obtain the frequency co c of the 
collective mode near the critical point as 



(18) 
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The constant £ can be determined by a variational 
method. According to Ref. |Q, an upper bound of 
the frequency of a collective mode is given by w^PP or = 
^/m 3 /m 1 /h, where m 1 = (0\[F, [H, F]]|0)/2 and m 3 = 
(0\[[F,H],[H, [H,F]]}\0)/2 are the energy- weighted mo- 
ment and the cubic-energy-weighted moment of the dy- 
namic structure factor, with F being an excitation op- 
erator. For the monopole mode the excitation operator 
is given by F = Y%Zi °( x i + Hi + Z D- Applying this 
formula to our Hamiltonian which is given in the square 
brackets of Eq. (|l|) , we obtain Eq. (|l^) with £ = 1/4 and 
u c replaced by o;" ppcr . Stringari Q has pointed out that 
for the case of repulsive interaction this method gives re- 
sults in excellent agreement with numerically obtained 
exact frequencies. We expect that the same is true for 
the case of attractive interaction and identify tL>" ppcr with 
lu c . Substituting Eq. @ into Eq. (|is|) and using Eq. (|), 
we obtain 



^ = 160* (l-^ 



(19) 



We thus find that as the number of condensate bosons Nq 
approaches its critical value N c , the collective frequency 
u! c vanishes as the one- fourth power of 1 — Nq/N c . 

Rate of Macroscopic Quantum Tunneling 

At sufficiently low temperature, the thermally acti- 
vated decay of the condensate over the barrier is negligi- 
ble, and we can expect MQT to provide a dominant decay 
mechanism. When the barrier is low enough for MQT 
to occur but still so high that the instanton approxima- 
tion is valid, the MQT rate T is given by T = Ae~ sB / h , 
where S B is the bounce exponent, that is, the value of 
the imaginary-time action S evaluated along the bounce 
trajectory r B (r). Using Eq. (jl4|) with £ = 1/4 we obtain 



S/h 



No 



dr[3r 2 + /(r)-/(r min )] , 



(20) 



where the imaginary time has been rescaled as lot — > r. 
The bounce trajectory r B (r) is the one that makes S 
extremal. ^From SS/5r B — 0, we obtain an equation of 
motion for the bounce trajectory which can be integrated 
to give (f B ) 2 = ±[/(r B ) - /(r min )]. Near the critical 
point, it is sufficient to expand the right-hand side of 
this equation in powers of r B — r mm and keep terms up 
to the third power: 



/(r B) _ /(r min ) 



(r 1 



„min \ 2 
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where r L is the left turning point of f(r) such that 
/(r L ) = f(r mm ). The bounce trajectory is then obtained 
as 



r B (r) = r" 



cosh 2 x //"(r min )/24r 



(22) 



Substituting Eqs. (Ell) and (E2) into Eq. (EOJ), we obtain 



S B /H = —N oy /6f"(r™ ia )(r min - r L ) 2 . (23) 
1 o 

Since r min - r L = 3/"(r min )//"'(r min ), we find from Eq. 
(0) and /"'(r min ) ~ 120/r c that r min -r L = 36r. Hence 



(24) 



Substituting Eqs. @ and (||) into Eq. @, we obtain 



(25) 



Thus if the condensate is formed and its decay is governed 
by MQT, the bounce exponent should be proportional to 
the five- fourths power of 1 — Nq/N c , in contrast to the 
unit power found in Ref. 

For a quadratic-plus-cubic potential, the prefactor A 
isgiven byA = w c (155 b /2tt?i)2 Substituting Eqs. 

( |l9| ) and (|25|) into this, we obtain 



- ~ 11.8aJ 



Ao 
A r 



(26) 



Thus as Nq approaches A c , the prefactor vanishes as the 
seven-eights power of 1 — No/N c . The crucial observation 
here is that the decrease in the prefactor as Ao — > A" c is 
outweighed by the much faster increase in the exponen- 
tial factor. Because of this rapid growth in the MQT 
rate, MQT can be a dominant decay mechanism of the 
condensate near the critical point for experiments by the 
Rice group || as we now show. 

Taking the experimental data of the second paper of 
Ref. ||, we have ui = (u) x u) y u) z )i ~ 908.4/sec, lo c ~ 
3231(1 - No/Nc) 1 / 4 , S B /h ~ 4.58AT (1 - N /N c ) 5 / 4 , and 
A ~ 10720Af Q 1/2 (l - Ao/A^ c ) 7 / 8 with N c = 1250. For 1 - 
Nq/N c = 10~ 2 , we obtain lu c ~ 1022/sec, A ~ 6707/sec, 
S B /h ~ 17.9, and T ~ 1.12 x 10~ 4 /sec. For this A 
MQT is negligible, but if Nq is a little bit closer to A^ c , 
e.g., for 1 — Nq/Nc — 5 x 10~ 3 , we obtain w c ~ 859/sec, 
A ~ 3666/sec, S B /h ~ 7.58, and T = 1.88/sec, and the 
MQT rate is therefore significant. 

It has sometimes been argued that the decay rate of 
the condensate due to MQT is much slower than that 
due to two-body dipolar and three-body collisions, and 
is therefore unlikely pjj. The above numerical evalu- 
ation, however, shows that near the critical point the 
MQT rate is at least comparable to the decay rate due 
to those inelastic collisions evaluated in Ref. |l7],[l£| be- 
cause near the critical point the MQT rate grows at an 
enormous rate as numerically illustrated above. Kagan 
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et dl. pointed out yet another interesting decay mecha- 
nism due to exchange interaction || . This mechanism be- 
comes important when the mean-field interaction energy 
per particle is larger than the single-particle energy-level 
spacing. For the experiments of Ref. 0, these energies 
are estimated to be InK and 7nK, respectively, so the 
condensate is not likely to decay via this mechanism. 

In conclusion, we have used a variational method and 
the instanton technique to find analytically the frequency 
of the collective mode and the MQT rate of a Bose con- 
densate with attractive interaction near the critical point 
of collapse. The maximum number of condensate bosons 
is found to be attained for the case of isotropic confining 
potential. By comparing MQT with other decay mech- 
anisms, we have argued that MQT can be a dominant 
decay mechanism of the condensate for Nq very close to 
N c . The obtained formulas contain no fitting parame- 
ters and can therefore be used as rather stringent tests 
for the existence of BEC. The fact that we have demon- 
strated the tunneling of the condensate to indefinitely 
low energy to occur does not necessarily imply that the 
physical system actually collapses, because we have not 
taken into account any effect of higher-order interactions. 
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